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The majoron provides an attractive dark matter candidate, directly associated to the 
mechanism responsible for spontaneous neutrino mass generation within the standard model 
SU(3) C <8> SU(2)l ® U(1)y framework. Here we update the cosmological and astrophysical con- 
straints on majoron dark matter which come from the Cosmic Microwave Background (CMB) and 
a variety of X- and gamma-Ray observations. 



PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 

I. INTRODUCTION 

It is by now established that only a small fraction (less 
than 20%) of the matter in the Universe is in the form 
of ordinary - i.e., baryonic - matter, while the rest is 
in the form of so-called "dark matter". The existence 
of dark matter is inferred by gravitational anomalies at 
very different scales, ranging from galactic scales to clus- 
ter scales, and all the way up to the cosmological scales. 
In particular, the recent 9-ycar data from the Wilkinson 
Microwave Anisotropy Probe (WMAP) [lj [2] have pro- 
vided even stronger support to the six-parameter ACDM 
model, although small-scale experiments like ACT [5] and 
SPT [3] hint to interesting, albeit discordant, deviations 
from this simple picture, like e.g. the presence of other 
relativistic degrees of freedom, in addition to the stan- 
dard model neutrinos, or deviations from ordinary grav- 
ity. The soon-to-be-released data from the Planck satel- 
lite will allow to measure the ACDM parameters with 
high precision, as well as to assess the significance of var- 
ious extension to the simplest model. 

If the ACDM model is certainly a phenomcnological 
success, nevertheless it is puzzling from the theoretical 
point of view in many aspects. On one side, the na- 
ture of both dark matter and dark energy, that together 
make up for more than 95% of the total energy budget 
of the Universe, is still unknown. On the other side, the 
theory of inflation, that explains the formation of the 
primeval seeds for density fluctuations from which galax- 
ies originate, is still waiting to be embedded in a more 
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fundamental theory. All these puzzles hint, and possibly 
have their solution in, some physics beyond the standard 
model (SM) of particle physics, or maybe in some modi- 
fication of general relativity. 

Although its precise nature is still unknown, there is 
no shortage of candidates for the role of dark matter. 
One of the most widely studied candidates to date is 
the supersymmetric neutralino. Recent results from the 
Large Hadron Collider (LHC), however, have greatly re- 
duced the available parameter space for supersymmetry, 
at least in its simplest minimal supergravity implementa- 
tions [5], reducing, perhaps, the appeal of supersymmet- 
ric dark matter candidates. Other possibilities include 
the axion, sterile neutrinos, Kaluza-Klein dark matter, 
and many others. 

The evidence for the existence of dark matter is very 
strong, but only limited to effects related to its grav- 
itational interaction. The search for non-gravitational 
evidence of the dark matter continues, in the form of 
direct and indirect detection experiments, and by look- 
ing for dark matter production in accelerators like the 
LHC. A precise underpinning of the dark matter and de- 
termination of its properties can only come through a 
combination of these approaches. 

If the dark matter has any connection to the world 
of SM particles, there will be astrophysical signals one 
can search for, in particular high energy photons from 
annihilating or decaying dark matter (see Ref. [5] for a 
recent review) . The most studied scenarios are the broad 
spectrum annihilation signals from neutralinos, but the 
real smoking gun is line emission (either directly from the 
decay /annihilation 7, 8 J or from internal bremsstrahlung 
0[in])j for which the spectral and spatial distribution is 
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not easily mimicked by astrophysical sources. The recent 



claim of line emission at E~ 



130 GeV in the Fermi 



data [111 1 1 2] has spurred a renewed interest in emission 
line searches at high energy. However, the origin of the 
possible signal at E~ f = 130 GeV is still unknown and 
requires further investigation. 



It has been long suggested [13] that the origin of 
dark matter could be related to the origin of neutrino 
masses [2]. In fact, the smallness of neutrino masses, 
as compared to the other SM particles, is puzzling in it- 
self. Most likely it is associated to the properties of the 
messenger states whose exchange is responsible for in- 
ducing them. This is the idea underlying the so-called 
seesaw mechanism [15H19j . whose details remain fairly 
elusive. Especially appealing is the possibility that neu- 
trino masses arise from the spontaneous violation of un- 
gauged lepton number [2"Ul |2"T] . The associated Nambu- 
Goldstone boson, the majoron, could acquire a mass from 
non-perturbative gravitational effects [23J [32] , and play 
the role of the dark matter particle. In Ref. [23] the vi- 
ability of the majoron as a dark matter particle was ex- 
plored using the WMAP 3-year data and in Ref. [25] the 
possible X-ray signature associated to majoron decay was 
investigated. A specific theoretical model implementing 
the seesaw mechanism and an A4 flavour symmetry was 
described in |21H. 



In this paper, we update our previous constraints 
in the light of the more recent cosmological and as- 
trophysical data. Regarding cosmology we take into 
account the most recent WMAP 9-year data [TJ [SJ. 
On the astrophysical front we include emission line 
searches on the entire range of photon energies between 
0.07 keV and 200 GeV from Chandra X-ray Observatory, 
X-ray Multi-Mirror Mission - Newton (XMM), High En- 
ergy Astronomy Observatory Program (HE AO), INTEr- 
national Gamma-Ray Astrophysics Laboratory (INTE- 
GRAL), Compton Gamma Ray Observatory (CGRO), 
and the Fermi Gamma-ray Space Telescope. 



The paper is organized as follows. In Sec [IT] we briefly 
recall the basics of majoron physics.. In Sees. |III| and [TV] 
we derive observational constraints on the majoron decay 
to neutrinos and photons, respectively, and we compare 
them to the predictions of a general seesaw model. Fi- 
nally, in Sec. [V]we draw our conclusions. 



II. SEESAW MAJORON PHYSICS 

The basic idea is that the lepton number symmetry 
of the standard SU(3) C ® SU(2) L ® U(1) Y model is pro- 
moted to a spontaneously broken symmetry [201 EI]- This 
requires the presence of a lepton-number-carrying com- 
plex scalar singlet coupling to the singlet neutrinos 
as follows, 



A<72/£ <T 2 2/£ + H.c. 



(1) 



This term provides the large mass term in the seesaw 
mass matrix 



V3W3 Y v v 2 
Y v T v 2 Y lVl 



(2) 



in the basis vl, f£ of "left" and "right" -handed neutrinos. 
The model is characterized by singlet, doublet and triplet 
Higgs scalars whose vacuum expectation values (vevs) are 
arranged to satisfy v± ^> v 2 3> v 3 obeying a simple vev 
seesaw relation of the type 



V 3 Vx ~ i> 2 



(3) 



The vev v\ drives lepton number violation and induces 
also a small but nonzero V3 while v 2 2 is fixed by the 
masses of the weak gauge bosons, the W and the Z. Note 
that the vev seesaw condition implies that the triplet vev 
V3 — > as the singlet vev v± — > 00. The three vevs deter- 
mine all entries in the seesaw neutrino mass matrix. Re- 
garding the Yukawa couplings, Y v is an arbitrary flavour 
matrix, while Y3 and Y\ are symmetric. The effective 
light neutrino mass obtained by perturbative diagonal- 
ization of Eq. ^ is of the form 

m v ~ Y 3 v 3 - Y V Y X - X Y V T — (4) 

Together with Eq. ^ this summarizes the essence of the 
seesaw mechanism. 

In order to identify which combination of Higgs fields 
gives the majoron, one may write the scalar potential ex- 
plicitly, minimize it, and determine the resulting scalar 
mass matrices. However, one can do this simply by ex- 
ploiting the invariance properties of the Higgs potential 
V [21]. The result is proportional to the combination 

J ex v 3 v 2 2 %(A )-2v 2 v 3 2 Z(<P ) + v 1 (v 2 2 +W)%(<j) (5) 

up to normalization (note 3 denotes imaginary parts). 
We remark the presence of the quartic lepton-number- 
conserving term 



$ t Ar 2 $*fj* + H.c. 



(6) 
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in the scalar potential. Here r 2 is the weak isospin Pauli 
matrix, and w 2 = v i = (c}> v 3 = (A). This term 
illustrates the need for mixing among neutral fields be- 
longing to all three Higgs multiplets in the expression 
for the majoron, Eq. ([5]). This implies that the majoron 
has an explicit coupling to two photons that leads to the 
possibility of indirect detection of majoron dark matter 
by searching for the corresponding high energy photons 



which we treat in Sec. IV 



We now turn to the form of the couplings of the ma- 
joron within the above seesaw scheme, characterized by 
spontaneous lepton number violation in the presence of 
singlet, doublet and triplet Higgs scalars. Again one can 
derive the form of the couplings of the majoron using 
only the symmetry properties, as described in Ref. [21] , 

~2^ V i 9ijV2Vj+H.C. 



£viik = 



(7) 



The result is a perturbative expansion for the majoron 
couplings 

TTli „ 



9ij 



vi 



(8) 



where the dots . . . denote higher order terms. One sees 
that, to first approximation, the majoron couples to the 
light mass-eigenstate neutrinos inversely proportional to 
the lepton number violation scale v\ = (a) and propor- 
tionally to their mass. With this we can compute the 
dark matter majoron decay rate to neutrinos as 



Jvv 



mj Ei(» 



(9) 



32tt 2vf 

Moreover, there is a sub-leading majoron decay mode to 
gammas. Within the general seesaw model this decay is 
induced at the loop level, resulting in [25] 



J 71 



64tt 3 



/ 



' N f 



2v 2 
1 V2V1 



2T. 



12m 2 f 



(10) 



where Nf, Qf, T£ and rrif denote respectively the color 
factor, electric charge, weak isospin and mass of the SM 
electrically charged fermions /. We note that this for- 
mula is an approximation valid for mj <C rrif, however 
we will always use the exact formula in the actual calcu- 
lations. 

The decay of the majoron dark matter to neutrinos 
provides the most essential and model-independent fea- 
ture of the majoron dark matter scenario, namely, it is a 
decaying dark matter model where the majoron decays 
mainly to neutrinos, a mode that is constrained from the 



CMB observations, as we discuss in Sec. Ill 



III. CMB CONSTRAINTS ON THE INVISIBLE 
DECAY J -> vu 

We start by deriving constraints on the majoron prop- 
erties from CMB anisotropy data. The majoron differs 
from most dark matter candidates in that it is unstable, 
since it must decay to neutrinos, as seen in Eq. Q, al- 
though it obviously has to be long-lived enough to be the 
dark matter. 

In order to investigate the observable effects of majoron 
decay on the CMB, the Boltzmann equation describing 
the phase-space evolution of DM particles must be mod- 
ified accordingly, as shown e.g. in Ref. [37], both at 
the background and at the perturbation level. The main 
effect of the late DM decay to invisible relativistic par- 
ticles is the increase of the late integrated Sachs- Wolfe 
(LISW), caused by the presence of an extra radiation 
component at small rcdshifts. This is reflected in an in- 
creased amount of power at the largest angular scales 
(i.e., small multipoles). Too large a decay rate would pro- 
duce too much radiation and too much large-scale power, 
and thus be at variance with observations. In Ref. [24] 
two of us have used this effect to constrain the majoron 
lifetime. However, we did not take properly into account 
the effect of majoron decay on the age of the Universe; 
this lead to an underestimate of the upper limit of the 
majoron decay rate. We have now corrected this. 

We use a modified version of CAMB [28], taking into 
account the finite lifetime of the majoron, to evolve the 
cosmological perturbations and compute the anisotropy 
spectrum of the CMB for given values of the cosmolog- 
ical parameters. In order to compute credible intervals 
and sample the posterior distributions for the parameters 
of the model, given some data, we use the Metropolis- 
Hastings algorithm as implemented in CosmoMC [29] 
(interfaced with our modified version of CAMB). We 
take uniform priors on all the parameters of the model, 
namely the present density parameters fl^h 2 and Q,^ m h 2 
of baryons and dark matter respectively, the Hubble con- 
stant Ho, the optical depth to recombination r rcc , the 
spectral index n s and amplitude A s (evaluated at the 
pivot scale ko = 0.002 Mpc -1 ) of the spectrum of pri- 
mordial scalar fluctuations, and the decay rate Tj^ v n of 
majoron to neutrinos. We marginalize over the ampli- 
tude of the contamination from the Sunyaev-Zel'dovich 
signal. 

We perform our analysis using the most recent WMAP 
9-year temperature and polarization data [HE]. In par- 
ticular, in temperature we include data up to ^ max = 
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FIG. 1. Two-dimensional WMAP-9 constraints on the ma- 
joron dark matter present density and decay rate to neutri- 
nos. The light (dark) shaded regions correspond to 68% (95%) 
credible intervals. 



1200. We use the latest (V5) version of their likelihood 
code, publicly avaiilable at the lambda websittQ 

For the purpose of the present analysis, we are mainly 
interested in the limits on the decay width of majoron to 
neutrinos, Tj^ vv . We get the following upper limit at 
95% C.L. 



Tj^ vv < 6.4 x lO^^s" 1 



(11) 



resulting in a lower limit to the majoron lifetime tj > 
50 Gyr, roughly four times the age of the Universe. We 
have also checked that this limit does not substantially 
change when we allow for the possibility of extra degrees 
of freedom at the time of recombination. 

In Fig. [I] we show 68% and 95% credible intervals in 
the (rj^„„, fldmh 2 ) plane. There is an evident anticor- 
relation that is explained by the fact, already discussed 
in Ref. [21], that the CMB is mainly sensitive to the 
amount of dark matter at early times (through the height 
of the first peak in the anisotropy spectrum) . For a given 
amount of dark matter in the early Universe, a larger de- 
cay rate would result in a smaller present amount of dark 
matter. 



FIG. 2. 3<r line emission constraints on the decay rate into 
two mono-energetic photons. Applies to all dark matter can- 
didates with this signature. The constraints are taken from 
yellow [30], orange [31] (conservatively rescaled by a factor 
of two due to mass estimate uncertainties as recommended in 
[32] ) , red [33] , grey [31] , purple [35] , blue [31] , cyan [33 , green 



IV. X- AND 7- RAY CONSTRAINTS ON THE 
PHOTON DECAY J -> 77 

One of the most interesting features of sponta- 
neous lepton number violation within the general 
SU(3) C <g> SU(2)l <8) U(l)y seesaw model is that the neu- 
trino decay mode in Eq. ^ is accompanied by a two- 
photon mode, Eq. (10 1, as a result of the Eq. This is 
constrained by a number of astrophysical observations. 
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A. Existing constraints 

In Fig. [2] we plot the emission line constraints over the 
wide range of photon energies of 0.07 keV to 200 GeV. 

The very soft X-ray emission is covered by Chandra 
Low Energy Transmission Grating (LETG) observations 
of NGC3227 (0.07-4.1 keV) [30] and a rocket borne light 
cryogenic spectrometer (0.25 — 1.1 keV) [3T] I39j. 

The 0.3 — 12 keV range is well covered with constraints 
from various objects observed with the Chandra and 
XMM X-ray telescopes [3lH3l l4tM8] In Fig. [2] we have 
chosen the strongest robust constraints from XMM obser- 
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vations of the Milky Way and M31 [35] and from Chandra 
observations of the Draco dwarf galaxy [34j^| 

The diffuse X-ray background observed with HEAO 
was searched for line emission by Ref. [35] over the 
range 3 — 48 keV. Line emission constraints have been 
derived from INTEGRAL SPI observations of the soft 7- 
ray background (20keV— 7MeV) (35]. For energies above 
those covered by INTEGRAL the constraints are two or- 
ders of magnitude worse as this range is only covered by 
a combination of the rather old COMPTEL and EGRET 
instruments onboard the CGRO. However, line emission 
constraints have been derived up to 100 GcV [37 . The 
most recent flagship for 7-ray searches is the Fermi 7-ray 
Space Telescope, for which line emission searches have 
been performed for the range of 7 — 200 GcV [38 . 

B. Future sensitivities 

The constraints in Fig. [2] can be improved by increas- 
ing the statistics or the spectral resolution. Increas- 
ing statistics (either by exposure time or by sensitiv- 
ity) improves the constraints as T^ cw = ^N°x/N™ W T° X , 
where 7V° x,now are the existing and new total of pho- 
tons per bin (assuming both source and background 
counts increase by the same amount). Increasing the 
spectral resolution improves the constraints directly as 
r „ew = E£™ HM /E™ WRM r™ and is consequently prefer- 
able but also technically more challenging. 



C. Model comparison 

We now compare the observational constraints ob- 
tained in the previous section to the predictions of differ- 
ent realizations of a general majoron seesaw model. In 
particular, we perform a random scan over the Yukawa 
matrices and vevs that characterize the seesaw mass ma- 
trix M. v in Eq. ([2]). For each point in the parameter 
space we evaluate the effective light neutrino mass ma- 
trix and the Majoron decay rate to neutrinos following 



2 Some analyses have claimed stronger constraints, but later found 
to be too optimistic. Ref. 1491 underestimated the flux by two 
orders of magnitude [401142] . According to Ref. [31T the mass was 
overestimated in Ref. 1501 leading to too restrictive constraints. 
The constraints in Ref. |51l might be too restrictive due to the 
choice of source profile 36 , and the spectral resolution appears 
overestimated in Ref. 1521 
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FIG. 3. The line emission constraints from Figure [2] (grey) 
compared to model predictions (colored dots), for different 
values of the triplet vev V3. 

Eqs. Q and ([9]), respectively. We then choose, among 
all possible realizations, those that are in agreement with 
current neutrino oscillation data [53] as well as with the 



neutrino bound in Eq. ( 11 ). Finally we compute the cor- 



responding decay rate to photons, as described in Sec. 

US 

We show the results of our scan in parameter space 
in Fig. [3j together with the constraints already shown 
in Fig. [2j It is clearly visible that the J — > 77 con- 
straints from line emission searches already begin to cut 
the parameter space for realistic models. This happens 
in particular for models with W3 larger than a few MeVs. 
However, models with lower values of the triplet vev pre- 
dict a photon flux that below the observational limits by 
at least 6 orders of magnitude. 



V. CONCLUSIONS 

We have updated previous constraints on the param- 
eters of the majoron dark matter model using the most 
recent CMB, X- and 7-ray observations. From the CMB, 
we have derived an upper limit on the rate of the invisible 
decay of the dark matter particle, namely, in the frame- 
work of the model under consideration, on the majoron 
decay to neutrinos. Translated in terms of the particle 
lifetime, this constraints the majoron lifetime to be larger 



6 



than 50 Gyrs. 



in parameter space can be reduced. 



The majoron also possesses a subleading decay mode 
to two photons, that can be constrained by astrophysical 
observations in the X and 7 regions. We have compared 
these limits to the theoretical predictions corresponding 
to different values for the parameters of the underlying 
particle physics model. We have found that the observa- 
tional constraints already exclude part of the parameter 
space for models in which the vev of the triplet W3 is 
larger than a few Me Vs. On the other hand, for smaller 
values of W3, the current limits need to be improved by 
at least 6 orders of magnitude before the allowed region 
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